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The quasi-elastic contribution of the nuclear inclusive electron scattering model de- 
veloped in A. Gil, J. Nieves, and E. Oset: Nucl. Phys. A 627 (1997) 543; is extended to 
the study of electroweak Charged Current (CC) induced nuclear reactions at intermediate 
energies of interest for future neutrino oscillation experiments. The model accounts for 
long range nuclear (RPA) correlations, Final State Interaction and Coulomb corrections. 
RPA correlations are shown to play a crucial role in the whole range of neutrino energies, 
up to 500 MeV, studied in this work. Predictions for inclusive muon capture for different 
nuclei, and for the reactions 12 C(^ M , /J-~)X and 12 C(v e , e~)X near threshold are also given. 
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1 Introduction 

There are two main reasons that motivate the study of neutrino scattering off 
nuclei. Firtsly the topic is of interest by itself because of the interest of knowing 
how the nuclear medium affects the nucleon axial current. The other main reason is 
that projected neutrino experiments require of a good control over systematic errors 
coming from the uncertainty in the neutrino-nucleus cross section. This second issue 
has motivated big efforts in this topic at low (few MeV) and intermediate nuclear 
excitations energies (a few hundred MeV). 

Any model aiming at describing the interaction of neutrinos with nuclei should 
be firstly tested against the existing data on the interaction of real and virtual 
photons with nuclei. The model developed in pp (inclusive electro-nuclear reac- 
tions) and |21 (inclusive photo-nuclear reactions) has been successfully compared 
with data at intermediate energies (nuclear excitation energies ranging from about 
100 MeV to 500 or 600 MeV). The building blocks of this model are: i) a gauge in- 
variant model for the interaction of real and virtual photons with nucleons, mesons 
and nucleon resonances with parameters determined from the vacuum data, and ii) 
a microscopic treatment of nuclear effects, including long and short range nuclear 
correlations 3 , Final State Interactions (FSI), explicit meson and A(1232) degrees 
of freedom, two and three nucleon absorption channels, etc. Nuclear effects are com- 
puted starting from a Local Fermi Gas picture of the nucleus, which is an accurate 
approximation to deal with inclusive processes which explore the whole nuclear 
volume, see [2] and 0]. The parameters of the model are completely fixed from 
previous hadron- nucleus studies: pionic atoms, elastic and inelastic pion- nucleus 
reactions, A— hypernuclei, etc. jo]. The gauge boson (real and virtual photons and 
charged bosons W^) coupling constants are also determined in the vacuum. Thus 
the model of £Q and [2] has no free parameters, and hence these results are predic- 
tions deduced from the nuclear microscopic framework developed in [3] and |H] . In 
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this talk, we show an extension of the nuclear inclusive quasi-elastic (QE) electron 
scattering model of pQ, including the axial CC current, to describe neutrino and 
antincutrino induced nuclear reactions in the QE region. We will not show here 
many details of the model, for a detailed discussion we refer the reader to 0. For 
an extension to neutral currents see [7] , where nucleon knock-out reactions are also 
estudied. 



2 INCLUSIVE CROSS SECTION 



We will expose here the general formalism focusing on the neutrino CC reaction. 
The generalization to antincutrino CC reactions or muon capture is straightforward. 
In the laboratory frame, the differential cross section for the process vi(k) + Az — * 




Fig. 1. Diagrammatic representation of the neutrino selfenergy in nuclear matter. 
l~(k') + X reads: 

lk ' IG \^ (!) 



<k!(£')ae; M 

with L and W the leptonic and hadronic tensors, respectively. On the other hand, 
the inclusive CC nuclear cross section is related to the imaginary part of the neu- 
trino self-energy (see Fig. ^) in the medium by: 



^ J d 3 rlm^(k; p(r)) (2) 



We get ImE„ by following the prescription of the Cutkosky's rules: in this case we 
cut with a vertical straight line (see Fig.Q the intermediate lepton state and those 
implied by the VF-boson polarization. Those states are placed on shell by taking 
the imaginary part of the propagator , self-energy, etc. We obtain for k° > 
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and thus, the hadronic tensor is basically an integral over the nuclear volume of 
the W— selfenergy (lJ^(q;p)) inside the nuclear medium. We can then take into 
account the in- medium effects (such as W— absorption by pairs of nucleons, pion 
production, delta resonances,...) by including the correspondent diagram in the 
shaded loop of Fig.^ Further details can be found in Ref. [JJ. 



3 QE CONTRIBUTION TO U^(q;p) 

The virtual W + can be absorbed by one nucleon (neutron) leading to the QE 
peak of the nuclear response function. Such a contribution corresponds to a lplh 
nuclear excitation. At sufficiently high energies other channels (such as pion pro- 
duction, two nucleon absorption, etc. ) open and should be taken into account. We 
consider a structure of the V — A type for the W + pn vertex, and use PC AC and in- 
variance under G-parity to relate the pseudoscalar form factor to the axial one and 
to discard a term of the form (p^ +7/^)75 in the axial sector, respectively. Invariance 
under time reversal guarantees that all form factors are real. Besides, and thanks to 
isospin symmetry, the vector form factors are related to the electromagnetic ones. 
We find 

f*rW) " " e W (r) - IPI) 



2M 2 J tf£(p)£(p+q) v * w ™ (4) 

6(|p + q| - k F (r))6(q° + E(p + q) - E(p))A^(p, q)\ p o =E{p) 

with the local Fermi momentum kp{r) = (3n 2 p(r)/2) 1 / 3 , M the nucleon mass, 
and E(p) = yj M 2 + p 2 . We will work on an non-symmetric nuclear matter with 
different Fermi sea levels for protons, k p F , than for neutrons, k F (equation above, 
but replacing p/2 by p p or p n , with p = p p + p n ). Finally, is the CC nucleon 
tensor 0. The d 3 p integrations above can be done analytically and all of them 
are determined by the imaginary part of isospin asymmetric Lindhard function, 
U(q, k F ,k F ). Explicit expressions can be found in UJ. 

We take into account polarization effects by substituting the particle-hole (lplh) 
response by an RPA response consisting in a series of ph and A-hole excitations 
as shown in Fig. We use an effective Landau-Migdal ph-ph interaction [H]: 
V = co {/o + /0T1T2 + go"" if 2 + 50C1C2T1T2}, where only the isovector terms 
contribute to CC processes. In the S = 1 = T channel (ctcttt operator) we use an 
interaction with explicit tt (longitudinal) and p (transverse) exchanges. The 

A(1232) degrees of freedom are also included in the S = 1 = T channel of the RPA 
response. This effective interaction is non-relativistic, and then for consistency we 
will neglect terms of order 0(p 2 /M 2 ) when summing up the RPA series. 

We ensure the correct energy balance, of both neutrino and antineutrino CC 
induced process in finite nuclei, by modifying the energy conserving 5 function in 
Eq. Q to account for the Q- value of the reaction. 

We also consider the effect of the Coulomb field of the nucleus acting on the 
ejected charged lepton. This is done by including the lepton self-energy Sc = 
2k'°Vc(r) in the intermediate lepton propagator of Fig.Q] 
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Fig. 2. Set of irreducible diagrams responsible for the polarization (RPA) effects in the 
lplh contribution to the W— selfenergy. 

Finally, we account for the FSI by using nucleon propagators properly dressed 
with a realistic self-energy in the medium, which depends explicitly on the energy 
and the momentum 9 . Thus, we rewrite the imaginary part of the Lindhard func- 
tion (ph propagator) in terms of particle and hole spectral functions S P! h(w, p; p)- 

4 RESULTS 

We present in Fig. [3] and Table ^ our theoretical predictions and a comparison 
of those to the experimental measurements of the inclusive 12 C (v^, [i~)X and 12 C 
{v e , e~)X reactions near threshold. Pauli blocking and the use of the correct energy 
balance improve the results, but only once RPA and Coulomb effects are included 
a good description of data is achieved. 



Table 1. Experimental and theoretical flux averaged 12 C(f M ,/i )X and 12 C(u e ,e )X 
cross sections in 10 -40 cm 2 units. We label our predictions as in Fig. |3] 





RPA 


Exp 10 and 11 




11.9 


LSND'95 
8.3 ±0.7 ±1.6 


LSND'97 
11.2 ±0.3 ±1.8 


LSND'02 
10.6 ±0.3 ±1.8 


a(u e , e~) 


0.14 


KARMEN 
0.15 ±0.01 ±0.01 


LSND 
0.15 ±0.01 ±0.01 


LAMPF 
0.141 ±0.023 



Given the success in describing the LSND measurement of the reaction 12 C 
{v^,^T)X near threshold, it seems natural to further test our model by studying 
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the closely related process of the inclusive muon capture in 12 C. Furthermore, 
and since there is abundant and accurate measurements on nuclear inclusive muon 
capture rates through the whole Periodic Table, we have also calculated muon 
capture widths for a few selected nuclei. Results are compiled in Table [5] Data 
are quite accurate, with precisions smaller than 1%, quite far from the theoretical 
uncertainties of any existing model. Medium polarization effects (RPA correlations) , 
once more, are essential to describe the data. Despite of the huge range of variation 
of the capture widths (note, r cxp varies from about 4xl0 4 s" 1 in 12 C to 1300 
xlO 4 s -1 in 208 Pb), the agreement to data is quite good for all studied nuclei, 
with discrepancies of about 15% at most. It is precisely for 12 C, where we find the 
greatest discrepancy with experiment. Nevertheless, our model provides one of the 
best existing combined description of the inclusive muon capture in 12 C and the 
LSND measurement of the reaction 12 C (v^, H~)X near threshold. 



Table 2. Experimental and theoretical total muon capture widths for different nuclei. 
Experimental data are taken from Ref. |12|. and when more than one measurement is 
quoted in .12', we use a weighted average: T/a 2 = X^R/°fj with 1/cr 2 = Y, 
We quote two different theoretical results: i) Pauli+Q obtained without including FSI 
effects and RPA correlations ; ii) the full calculation, including all nuclear effects with 
the exception of FSI, and denoted as RPA. In the last column we show the relative 
discrepancies existing between the theoretical predictions given in the third column and 

data. 





Pauli+Q [10 4 s _1 ] 


RPA [10 4 s- 1 ] 


Exp [lCs^ 1 ] 






L2 C 


5.42 


3.21 


3.78 ± 0.03 




0.15 


16 


17.56 


10.41 


10.24 ±0.06 




-0.02 


18 


11.94 


7.77 


8.80 ±0.15 




0.12 


23 Na 


58.38 


35.03 


37.73 ±0.14 




0.07 


40 Ca 


465.5 


257.9 


252.5 ± 0.6 




-0.02 


44 Ca 


318 


189 


179 ±4 




-0.06 


75 As 


1148 


679 


609±4 




-0.11 


112 Cd 


1825 


1078 


1061±9 




-0.02 


208 p b 


1939 


1310 


1311±8 




0.00 



At intermediate energies the predictions of our model become reliable not only 
for integrated cross sections, but also for differential cross sections. We present re- 
sults for incoming neutrino energies within the interval 150-400 (250-500) MeV for 
electron (muon) species. In Figs.0]and|3 FSI effects on differential cross section are 
shown. As expected, FSI provides a broadening and a significant reduction of the 
strength of the QE peak. Nevertheless the integrated cross section is only slightly 
modified. In Table|3]we compile electron neutrino and antineutrino inclusive QE in- 
tegrated cross sections from oxygen. Though FSI change importantly the differential 
cross sections, it plays a minor role when one considers total cross sections. When 
medium polarization effects are not considered, FSI provides significant reductions 
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E v -E(-Q [MeV] E v - E( - Q [MeV] 



Fig. 3. Theoretical predictions for the 12 C (y ll ,n~)X and the 12 C (v e ,e~)X reactions 
near threshold. In addition to the full calculation (denoted as RPA), with all nuclear 
effects with the exception of the FSI ones, we show results obtained without including 
RPA, FSI and Coulomb corrections (denoted as Pauli), and also results (denoted as LDT) 
obtained by multiplying the free space cross section by the number of neutrons of 12 C. 



(15-30%) of the cross sections. However, when RPA corrections are included the 
reductions becomes more moderate, always smaller than 7%, and even there exist 
some cases where FSI enhances the cross sections. This can be easily understood by 
looking at Fig. [5] There, we see that FSI increases the cross section for high energy 
transfer. But for nuclear excitation energies higher than those around the QE peak, 
the RPA corrections are certainly less important than in the peak region. Hence, 
the RPA suppression of the FSI distribution is significantly smaller than the RPA 
reduction of the distribution determined by the ordinary Lindhard function. 
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Fig. 4. Inclusive QE double differential cross sections in oxygen as a function of the 
transferred energy, for two values (310 and 450 MeV) of the transferred momentum. The 
incoming neutrino (antineutrino) energy is 400 MeV. We show results non-relativistic 
nucleon kinematics and with (FSI) and without (NOREL) FSI effects. For all cases, we 
show the effect of taking into account RPA correlations and Coulomb corrections (lower 

lines at the peak). 




Fig. 5. Muon neutrino inclusive QE differential cross sections in oxygen as a function of 
the transferred energy. The incoming neutrino energy is 375 MeV. The notation for the 
theoretical predictions is the same as in Fig. 2] 
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Table 3. Electron neutrino (left) and antineutrino (right) inclusive QE integrated cross sec- 
tions from oxygen. We present results for relativistic ('REL') and non-relativistic nucleon 
kinematics. In this latter case, we present results with ('FSF) and without ('NOREL') 
FSI effects. Results, denoted as 'RPA' and 'Pauli', have been obtained with and without 
including RPA correlations and Coulomb corrections, respectively. 



E v [MeV] 




a(">0(is e ,e-)X) 10- 




a ( ib O(/> e ,e+)X) KT 


- 4 w 






REL 


NOREL 


FSI 


REL 


NOREL 


FSI 


400 


Pauli 


389.4 


416.6 


352.5 


130.0 


139.1 


121.0 




RPA 


294.7 


322.6 


303.6 


91.9 


101.9 


104.8 


310 


Pauli 
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240.6 
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87.2 
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209.0 


195.2 
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72.4 


73.0 


220 


Pauli 


149.5 
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51.0 




RPA 
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